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Abstract 

Cooling of neutron stars (NSs) with superfluid cores is 
simulated taking into account neutrino emission pro- 
duced by Cooper pairing of nucleons. Two cooling 
regimes of NSs composed of matter with the same mod- 
erately stiff equation of state are studied: the stan- 
dard cooling (for a 1.30 Mq star, as an example) and 
the cooling enhanced by the direct Urea process (for 
a 1.48 Mq star). The critical temperatures of neutron 
and proton superfluidities, T cn and T cp , are assumed 
to be constant over the NS core, and treated as free 
parameters. They are constrained using the surface 
temperatures T s of isolated NSs (EX J0822-43, PSE 
1055-52, IE 1207-52, Vela, Geminga, PSE 0656+14, 
EX J0002+62), obtained by interpretation of observed 
thermal radiation either with black body spectrum or 
with hydrogen atmosphere models. 

The temperatures T s given by the "atmospheric" in- 
terpretation for the last five objects can be explained 
by cooling of the same NS, with the same T cn and 
T cp for all objects. The confidence regions of T cn and 
T cp are rather narrow and depend on cooling regime. 
Standard cooling requires moderately strong neutron 
and proton superfluidities, while enhanced cooling re- 
quires moderate neutron and strong proton superflu- 
idities. Adopting interpretation of observations with 
the black body spectrum, we can propose similar ex- 
planation only for three last objects. For the both, 
standard and enhanced, cooling regimes the confidence 
regions of T cni T cp appear to be wider than in the case 
of the "atmospheric" values of T s and consist of two 
separate subregions. None of the models requires si- 
multaneously strong neutron and proton superfluidities 
(T cn , T cp > 3 • 10 9 K), which is an argument against 
very soft equations of state in the NS cores. 



1 Introduction 

According to current theories, neutron stars (NSs) are 
born hot, and then gradually cool. The first 10 5 -10 6 yr 
cooling proceeds via neutrino emission from the NS 
core, but later it goes via photon surface emission. 
Cooling depends on NS internal structure. Thus, com- 
parison of cooling theory with observational data on 
NS surface temperatures enables one, in principle, to 
explore NS structure, particularly, properties of super- 
dense matter in the NS cores. Among various prop- 
erties of this matter, the critical temperatures of neu- 
tron and proton superfluidities, T cn and T cp , are of 
special interest. These quantities are sensitive to pa- 
rameters of strong nucleon-nucleon interaction in su- 
perdense matter closely related to the main "mystery" 
of NSs, the equation of state in their cores (soft, moder- 
ate or stiff?). The softer equation of state, the weaker 
nucleon-nucleon repulsion at small separations, and 
the stronger is nucleoli superfluidity (i.e., higher T cn 
and T cp ). The superfluidity affects the NS heat capac- 
ity and neutrino emission, and thus the stellar cool- 
ing. The aim of the present article is to continue our 
studies (e.g., Levenfish & Yakovlev 1996 and references 
therein) of cooling of NSs with superfluid cores. As in 
our previous articles, we assume, for simplicity, that 
the NS cores consist of neutrons (n), with admixture 
of protons (p) and electrons (e) . 

Superfluidity was introduced into the cooling theory 
by Tsuruta et al. (1972). However, it has not been re- 
garded as a powerful cooling regulator for a long time, 
and it has been taken into account in a simplified man- 
ner (as discussed, for instance, by Levenfish & Yakovlev 
1996). Page and Applegate (1992) were the first who 
emphasized the importance of superfluidity in cooling 
of massive NSs, where enhanced neutrino reactions of 
the direct Urea process are allowed: 
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The authors showed that the surface temperature of 
such a star was almost solely determined by the val- 
ues of T cn or T cp in its core (they assumed the super- 
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fluidity of one nucleoli species) in a wide interval of 
stellar ages, 10 3 - 10 5 yr. This means that NSs are 
good "thermometers" of nucleon superfluidity. More- 
over, the superfluidity appears to be a powerful regu- 
lator of standard cooling of less massive NSs (see, e.g., 
Page 1994, Levenfish & Yakovlev 1996) determined by 
much weaker neutrino reactions of the modified Urca- 
processes 

n + N -> p + N + e + D e , p+e + N ->• n+ N+v e , (2) 

and the nucleon-nucleon scattering 

N + N ^N + N + v + D, (3) 

where N is a nucleon (n or p). 

NS cooling is also greatly affected by neutrino emis- 
sion due to direct interband transitions of nucleons 
whose dispersion relation contains a superfluid gap 
(Schaab et al. 1997, Page 1998, Yakovlev et al. 1998): 

N -> N + v + v. (4) 

This process can be treated as neutrino generation due 
to Cooper pairing of (quasi) nucleons. It has been in- 
cluded into the cooling theory quite recently. Thus, 
previous cooling calculations have to be reconsidered. 

Theoretical conclusion on the importance of super- 
fluidity in NS cooling was made at about the same time 
when new observational data on NS surface tempera- 
tures appeared. Below we present new cooling calcu- 
lations (Sect. 2) and compare them with observations 
(Sect. 3). 

2 Cooling calculations 
2.1 Cooling models 

Our cooling simulations have been described by Lev- 
enfish & Yakovlev (1996). Here we point out some 
details. We use the models of NSs whose cores are 
composed of matter with moderately stiff equation of 
state proposed by Prakash et al. (1988) (the version 
with the compression modulus Kq = 180 MeV and 
with the symmetry energy Sy suggested by Page & 
Applegate 1992). The maximum NS mass, for a given 
equation of state, is 1.73 Mq. In order to study the 
enhanced and standard cooling we consider NS models 
with two masses M. In the first case M — 1.48 Mq, 
the NS radius is R — 11.44 km, and the central density 
/9 C = 1.376-10 15 g/cm 3 , while in the second case M = 
1.30 Mq, R = 11.87 km, and p c = 1.07-10 15 g/cm 3 . The 
adopted equation of state allows the direct Urea pro- 
cess to operate at densities p> p cr = 1.30 • 10 15 g/cm 3 . 
Thus the 1.48 Mq NS suffers the enhanced cooling: the 
powerful direct Urea process (pi) is allowed in a small 



central kernel of radius 2.32 km and mass 0.035 Mq, in 
addition to all the neutrino processes (§)-(Q) operat- 
ing in the entire core. On the other hand, the critical 
density p cr is not reached in the 1.30 Mq NS, and the 
star has the standard neutrino luminosity determined 
by the processes (||)-(||). Notice that while composing 
the equation of state in our previous cooling calcula- 
tions the parameter hq (saturation baryon number den- 
sity) has been set equal to 0.165 fm~ 3 . In the present 
simulations, we set hq — 0.16 fm~ 3 . That is why the 
NS mass for the enhanced cooling models is slightly 
different from that adopted earlier (1.44 Mq). 

Numerous microscopic calculations of Cooper pair- 
ing in the NS cores (which use different models of 
nucleon-nucleon interaction and manybody theories; 
see, e.g., Takatsuka & Tamagaki 1997 and references 
therein) predict a very large scatter of superfluid crit- 
ical temperatures of neutrons and protons, T cn and 
T cp (~ 10 7 -10 10 K), and their different density depen- 
dence. Therefore we will make a simplified assumption 
that the critical temperatures T cn and T cp are constant 
within the stellar core and can be treated as free pa- 
rameters. This is the main simplification of the present 
article. Let us assume that the proton superfluidity is 
produced by the singlet-state ( 1 So) pairing, while the 
neutron superfluidity is due to the triplet-state ( 3 P2) 
pairing with zero projection of total Cooper-pair mo- 
mentum onto quantization axis. We will study cooling 
for different T cn and T cp , and determine those T cn and 
T cp which are in better agreement with observations. 

The heat capacity and neutrino luminosity in the 
reactions (@)-(§) in the superfluid NS cores have been 
calculated as prescribed by Levenfish & Yakovlev (1996), 
while the luminosity in the reaction (^) has been taken 
from Yakovelv et al. (1998). We have also incorpo- 
rated the neutrino luminosity of the NS crust due to 
bremsstrahlung emission by electrons which scatter off 
atomic nuclei (using an approximate formula proposed 
by Maxwell 1979). The NS heat capacity has been de- 
termined as a sum of the capacities of n, p and e in 
the stellar core; we have neglected the heat capacity of 
the crust due to the smallness of the crust mass in our 
NS models. While calculating the neutrino luminosity 
and heat capacity the effective masses of neutrons and 
protons in stellar matter, to* and m*, have been set 
equal to 0.7 of their bare masses. 

As in our previous articles, we have used the NS 
cooling code based on the approximation of isothermal 
interior. The effects of General Relativity on NS struc- 
ture and cooling have been incorporated explicitly. The 
isothermal approximation is justified for a NS of age 
t > (10 - 10 3 ) yr, when the internal thermal relaxation 
is over, and only a narrow layer of the outer NS crust is 
non-isothermal. In the isothermal region the quantity 
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Figure 1: Effective surface temperature T s °° (red- 
shifted) versus age for standard cooling of a NS 
(1.30 Mq) with the superfluid core. Curves are la- 
beled by the values of lg T C n and lgT C p- Solid and 
dashed lines are calculated including and neglecting 
neutrino emission produced by nucleon pairing, respec- 
tively. The solid and dashed curves for \gT cn = 7.8 and 
lg T cp — 7.8 coincide. Dots show cooling of a nonsuper- 
fluid star. 



Ti(t) = T(r,t) ■ exp[$(r)] (which can be treated as the 
internal temperature corrected due to the local gravi- 
tational redshift) is constant at any moment of time t. 
Here, T(r,t) is the real local temperature, <i>(r) is the 
dimensionless gravitational potential, and r is a radial 
coordinate. 

The photon NS luminosity depends on the effective 
local surface temperature T s . The relationship between 
T s and the internal temperature Tj is determined by 
thermal insulation of the outer stellar envelope. We 
have ignored the effect of magnetic fields on NS cool- 
ing, and have used the relationship T s (Ti), obtained 
recently by Potekhin et al. (1997) for B = 0. We 
assume that the NS may possess a thin surface layer 
(of mass < lO _13 Af ) of hydrogen or helium. This 
layer cannot affect thermal insulation and cooling but 
can change spectrum of thermal emission. The effect 
of dipole surface magnetic fields i?<5-10 12 Gon 
NS cooling is not large (Shibanov & Yakovlev 1996). 
Therefore our calculations can be applied, at least qual- 



Figure 2: Enhanced cooling of a 1.48 M NS. No- 
tations are the same as in Fig. 1. The solid and 
dashed lines for lgT c „ = 8.3, lgT cp = 8.6 and for 
lgT c „ = 10.0, lgTcp = 8.4 coincide. Dots show cooling 
of a nonsuperfluid star. 

itatively, for the NSs with such magnetic fields. In 
these cases, T s means the mean effective surface tem- 
perature which determines the total photon NS lumi- 
nosity L 7 = 47rcri? 2 T s 4 ( non-redshifted) , where a is the 
Stefan-Boltzmann constant. 

2.2 Effect of neutrino emission produced 
by Cooper pairing of nucleons 

The main difference of the present simulations from our 
previous ones is in including the neutrino emission (Q) 
due to Cooper pairing of nucleons. The process was 
suggested by Flowers et al. (1976) and independently 
by Voskresensky & Senatorov (1987), but it has been 
included into the cooling simulations only recently (see 
above). The reaction rates have been recalculated and 
improved by Yakovlev et al. (1998). As shown by the 
latter authors, neutrino emission produced by proton 
pairing is strongly suppressed due to numerical small- 
ness of weak vector neutral currents of protons. The 
importance of the process for standard and enhanced 
NS cooling is illustrated in Figs. |l| and ||. 

It has been widely recognized that the superfluid- 
ity onset at the neutrino cooling stage slows down the 
standard cooling. As seen from Fig. [j], superfluidity, on 
the contrary, can accelerate the cooling. The accelera- 
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Figure 3: Lines of values of T cn and T cp which corre- 
spond to certain internal temperatures Tj (the values 
of lg Ti are given near the curves) or surface tempera- 
tures T s °° (given in Table of a NS with enhanced neu- 
trino luminosity (1.48 Mq) and Geminga's age (3.4- 10 5 
yr). The darker the background the cooler is the star. 
The region of joint neutron and proton superfluidity 
(in the center and the upper right corner) is enclosed 
by dashed lines. A small region where the superfluidity 
has not appeared at given age (the left lower corner) is 
separated by the thick solid line. 

tion can be so large that by the end of the neutrino era 
the surface temperature falls below the value it would 
have for the enhanced cooling. This conclusion does 
not apply to the cases of purely proton superfluidity 
(proton pairing produces weak neutrino emission), and 
strong purely neutron superfluidity (neutrino emission 
due to neutron pairing dominates over the standard 
neutrino processes @ and (|^) only for T < 10 9 K, see 
Yakovlev et al. 1998). 

Neutrino emission due to pairing of neutrons affects 
also enhanced cooling of NSs with superffuid protons 
at T cp ^ T cn (Fig. ||). In these cases, the proton su- 
perfluidity appears in the NS much earlier than the 
neutron one, and it has enough time to reduce the 
direct Urea process by the beginning of the neutron 
pairing. Therefore a splash of neutrino emission due 
to the neutron pairing appears to be quite substantial; 
the cooling does not slow down, as it would be without 
Cooper neutrinos, but strongly speeds up. 



Figure 4: Same as in Fig. | but for a NS of the 
PSR0656+ 14 pulsar age (10 5 yr). 

2.3 Results 

We have calculated about 2000 cooling curves which 
give the dependence of the effective NS surface tem- 
perature T s °° = T s y/l — Rg/R (as detected by a dis- 
tant observer, R g being the gravitational radius) on 
age t. For M = 1.30 M Q and M = 1.48 Af© we have 
T™/T a = 0.822 and 0.786, respectively. The critical 
temperatures of neutrons T cn and protons T cp in the 
NS cores have been varied in a wide interval from 10 6 
to 10 10 K, compatible with the scatter predicted by 
microscopic theories (see Sect. 2.1). 

Cooling curves for some selected values of T cn and 
T cp are presented in Figs. [I] and ||. However it is incon- 
venient to analyse the results in the form of the cool- 
ing curves. It is more convenient to plot those values 
of (T cn , T cp ), which lead to a selected surface tem- 
perature T s °° of the NS at given age t. Figures ||-|| 
are arranged in this manner. For example, we have 
chosen the values of t which correspond to the ages of 
the Geminga pulsar (3.4 • 10 5 yr, Figs. || and ||) and 
PSR0656 + 14 (10 5 yr, Figs. | and |). Observational 
data on these and other cooling NSs are given in Sect. 3. 

In order to simplify understanding of the figures let 
us outline the main properties of the heat capacity and 
neutrino emission from superfluid NSs (see Levenfish 
& Yakovlev 1996, Yakovlev et al. 1998, and references 
therein). The heat capacities of p and n constitute, re- 
spectively, ~ 1/4 and ~ 3/4 of the total heat capacity 
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Figure 5: Isotherms of the internal temperature Tj 
(or the surface temperature , Table [l]) of a NS of 
the Geminga's age, as in Fig. but for the standard 
cooling (1.30 M Q ). 

Ctot- Thus a strong proton superfluidity reduces Ctot 
by ~ 25%, while a strong neutron superfluidity reduces 
it by about 4 times. If the ratio of the inner temper- 
ature Tj to the critical nucleon temperature T c n in a 
cooling NS falls down from 0.3 to 0.1, the heat capac- 
ity of nucleons N = n or p decreases more than by 
three orders of magnitude, and becomes much lower 
than the heat capacity of electrons. Further reduction 
of the heat capacity of nucleons N does not affect the 
total heat capacity. The appearance of a weak superflu- 
idity of nucleons N almost doubles their heat capacity 
due to the latent heat release at the phase transition. 
The heat capacity of nucleons N remains higher than 
the heat capacity of normal nucleons when Ti/T c ^ de- 
creases from 1 to 0.3, i.e., as long as the difference of 
temperature logarithms is \gT C N — lgT; < 0.5. 

The effects of superfluidities of n and p on the neu- 
trino reactions are also different. For the reactions ([[])- 
(||) the difference is small. If we neglected neutrino 
emission due to Cooper pairing of nucleons the asym- 
metry of Figs. H"!! with respect to the inversion of axes 
T cn ^1 T cp would be mainly explained by different con- 
tributions of n and p into the heat capacity. Cooper- 
pairing neutrinos noticeably amplify the asymmetry: 
the neutrino emission due to pairing of p is weak, while 
the emission due to pairing of n dominates over the 



Figure 6: Same as in Fig. g, but for a NS of the 
PSR0656+ 14 pulsar age. 

standard neutrino energy losses at Tj < 10 9 K < T cn 
and over the direct Urea process at Ti < T cn <C T cp . 
The main energy release in the process occurs while 
Ti/T C N decreases from 0.96 to ~ 0.2 (when the differ- 
ence of temperature logarithms is \gT C N — lg < 0.7). 

Finally let us mention that for the enhanced cool- 
ing the transition from the neutrino to photon cool- 
ing stages occurs somewhat later than for the stan- 
dard cooling. In a nonsuperfluid NS which undergoes 
the standard cooling the neutrino era lasts t v ~ (1-3) 
•10 5 yr. A strong neutron superfluidity decreases t„ 
by a factor of 2-3, while a strong proton superfluidity 
decreases t v by about (20-30)%. If the both superflu- 
idities, n and p, are strong at once t v ~ 10 yr. 

Enhanced cooling 

Figures || and || illustrate enhanced cooling of NSs of 
ages of the Geminga and PSR0656 + 14 pulsars, re- 
spectively (see Table || below) . 

Figure || shows isotherms of the internal tempera- 
ture Ti of a NS of the Geminga's age versus T cn and 
T cp . Since Ti is strictly related to the surface temper- 
ature these lines are also isotherms of T s °° (Table Q). 
Dashes show auxiliary lines which enclose the region of 
the joint superfluidity of nucleons. The region to the 
left of the upper dashed line corresponds to superflu- 
idity of protons and normal neutrons, while the region 
to the right of the lower line corresponds to superfluid- 
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Table 1: Relationships between T t and Tf 3 [K] for Figs. |-|. 



M = 1.48 Mg 


IgTf 


5.02 5.28 5.54 5.79 6.03 6.25 6.36 6.40 6.46 6.56 6.58 6.60 6.62 
4.50 4.60 4.70 4.80 4.90 5.00 5.05 5.07 5.10 5.15 5.16 5.17 5.18 


IgT~ 


6.66 6.69 6.83 6.85 7.03 7.22 7.23 7.40 7.47 7.58 7.76 8.02 8.17 
5.20 5.22 5.29 5.30 5.40 5.50 5.51 5.60 5.64 5.70 5.80 5.95 6.03 



M = 1.30 Mq 


lgTf 


5.70 5.75 5.85 5.90 6.95 6.03 6.26 6.67 7.23 7.24 7.30 7.32 
4.77 4.79 4.83 4.85 4.87 4.90 5.00 5.20 5.49 5.50 5.53 5.54 


Wi 

lgT- 


7.37 7.39 7.43 7.50 7.55 7.61 7.64 7.68 7.70 7.75 7.77 7.82 
5.57 5.58 5.60 5.64 5.67 5.70 5.72 5.74 5.75 5.78 5.79 5.82 


lgTf 


7.84 7.86 7.88 7.91 7.93 7.98 8.00 8.05 8.07 8.13 8.14 8.18 
5.83 5.84 5.85 5.87 5.88 5.91 5.92 5.95 5.96 5.99 6.00 6.02 



ity of neutrons and normal protons. The dashed lines 
intersect at the isotherm of temperature Tj = 10 6,4 K, 
to which a nonsuperfluid star would cool down at mo- 
ment t. This isotherm, plotted by the thick line, en- 
closes the region where the nucleon superfluidity has 
not appeared by moment t and does not affect the 
cooling. Notice that owing to General Relativity effects 
isotherms Tj correspond to somewhat higher local tem- 
peratures of matter (see above). Thus, the isotherm 
T.- L = 10 6 4 K is associated with the values of T c m which 
are slightly higher than Ti. 

First of all we discuss the behaviour of isotherms 
to the right of the lower auxiliary line. The vertical 
parts of isotherms reveal that NS cooling is governed 
by the only superfluidity of n. With increasing T cn at 
lglcn <1 7.7, the core temperature Tj grows up which 
can be explained like this. The higher T cn , the ear- 
lier the neutron superfluidity appears; accordingly the 
powerful direct Urea process is suppressed earlier, cool- 
ing delay is longer, and the NS is hotter at given t. For 
low T cn , the NS cools down in a nonsuperfluid state 
the main part of its history, and is sufficiently cold at 
age t. The neutrino luminosity of such a star is rather 
weak and becomes comparable to the photon luminos- 
ity. Additional neutrino energy losses due to neutron 
pairing at lgT c „ > 6.6 slightly accelerate the cooling at 
6.6 < lgT c „ < 6.6 + 0.7. However at lgT c „ < 6.6 + 0.5 
the acceleration is compensated by the latent heat re- 
lease produced by the jump of the neutron heat capac- 
ity. The values 0.7 and 0.5 have been explained in the 



beginning of this section. 

Vertical parts of isotherms intersect the auxiliary 
line at temperatures T cpj at which the proton su- 
perfluidity appears. With further increase of T cp , a 
newly born proton superfluidity slightly decreases the 
neutrino luminosity which is already reduced by the 
strong neutron superfluidity The main effect will be 
caused by the heat capacity of protons which jumps up. 
The latent heat release will increase Tf. in a strip of 
width ^0.5 (in lgX" C p) above the lower dashed line the 
isotherms shift to the left. With further growth of T cp 
the heat capacity is strongly reduced, the NS becomes 
cooler and the isotherms shift to the right. 

In a similar fashion, exchanging proton and neutron 
superfluidities, we can explain horizontal parts of the 
isotherms to the left of the upper auxiliary line. In this 
case NS cooling in the range 6.6 < \gT cp < 6.6 + 0.7 
does not appear because of the weakness of the neu- 
trino emission produced by proton pairing. Strong pro- 
ton superfluidity suppresses the heat capacity weaker 
than neutron superfluidity with the same critical tem- 
perature. Therefore, at high T cp the star is warmer 
than for the same T cn at vertical parts of isotherms in 
the lower right part of the figure. 

The horizontal parts of the isotherms intersect the 
upper auxiliary line at temperatures T cn which switch 
on the neutron superfluidity. This superfluidity in- 
duces the latent heat release and associated slight NS 
heating (dips on isotherms to the right of the auxiliary 
line). With increasing T cp the NS heating produced 
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by neutron superfluidity decreases and disappears (the 
dips vanish). The effect is caused by the neutrino emis- 
sion due to neutron pairing which is more pronounced 
if the direct Urea process is strongly suppressed. At 
high T cp the neutron pairing becomes the main neu- 
trino emission mechanism. Since the neutrino emission 
due to neutron pairing is more efficient than the latent 
heat release the appearance of the neutron superfluid- 
ity with the growth of T cn does not delay the cooling 
but, on the contrary, accelerates it. A noticeable NS 
cooling via Cooper-pairing neutrinos, takes place in a 
strip of width 0.7 (in lgT cra ) to the right of the upper 
dashed line. Lower and to the right of this strip the 
cooling is associated with the reduction of the neutron 
heat capacity. 

Now consider the region of T cn , T cp between the 
auxiliary lines in more detail. Some increase of Tj at 
lgT c „ < 6.6 + 0.5 and lgT cp < 6.6 + 0.5 is caused by 
the latent heat release at onset of weak superfluidity of 
n and p. Further growth of T cn ~ T cp induces initially 
rapid decrease and then weak increase of Tj. The de- 
crease is explained by exponential reduction of the heat 
capacity by the joint superfluidity of n and p while the 
weak increase is associated with reduction of the direct 
Urea process at early cooling stages. 

Figure || is analogous to Fig. ||[ but corresponds to a 
younger NS of age t = 10 5 yr. Comparing Figs. || and 
[| we see that one should have T cp <C T cn or T cn <C T cp 
to support high surface temperature T s °° (or Ti) for a 
longer time. 

Standard cooling 

Standard cooling of a 1.30 Mq NS is illustrated in 
Figs. H and || for the stars of Geminga and PSR 
0656+ f 4 ages, respectively. Isotherms are qualitatively 
different from those for enhanced cooling (cf. with Figs. 
|^ and |J) : even an approximate symmetry of the neu- 
tron and proton superfluidities is absent. The asymme- 
try is attributed to the weakness of the standard neu- 
trino energy losses. First, in the absence of such a pow- 
erful cooling regulator as the direct Urea process the 
difference of the heat capacities of n and p (see above) 
is more pronounced. Second, Cooper-pairing neutrino 
emission becomes more important on the background 
of the weaker neutrino emission produced by other neu- 
trino reactions; Cooper-pairing emission is asymmetric 
itself being more efficient for neutrons than for protons. 

If the superfluidity is absent, a NS (1.30 M ) enters 
the photon cooling stage at t v ~ 1.6 • 10 5 yr. Thus the 
PSR 0656 + 14 pulsar appears at the transition stage 
and Geminga is at the photon stage. The neutrino 
luminosity is already weak and superfluidity affects 
the cooling mainly either through the heat capacity 
or through the Cooper-pairing neutrino emission. 



Consider, for instance, Fig. |5j. Two dashed auxil- 
iary lines enclose the domain of joint superfluidity of 
nucleons. To the left of the upper line protons are su- 
pcrfluid and neutrons not, while below the lower line 
only neutrons are superfluid. The lines intersect at the 
isotherm of the temperature lgTi ks 7.82, which a non- 
superfluid NS would have at age t. The superfluidity 
with lgT c „ < 7.9 and lgT cp < 7.9 has not appeared by 
moment t and does not affect the cooling. 

Horizontal parts of isotherms to the left of the up- 
per dashed line show that cooling is regulated by the 
proton superfluidity alone. At lg'T cp > 7.9 K the su- 
perfluidity appears just before given moment t and is 
weak. In the range 7.9 < lgT cp < 7.9 + 0.5 it initiates 
the latent heat release and a weak increase of Ti. At 
higher lgT cp > 8.4 the proton heat capacity is reduced, 
and the total NS heat capacity decreases by ~ 25%. 
On the other hand, for high T cp , the proton superfluid- 
ity onset is shifted to the neutrino cooling stage. The 
cooling delay produced by suppression of the neutrino 
luminosity in the neutrino era is somewhat stronger 
than the cooling acceleration produced by the effect of 
superfluidity on the heat capacity. Thus Tj continues 
its growth with increasing T cp . 

Horizontal parts of isotherms end at temperatures 
T cn ~ 10 8 K, at which the neutron superfluidity is 
switched on. The neutron pairing induces a splash of 
neutrino emission and cooling acceleration in a strip of 
the width ~0.7 (in lgT c „) to the right of the auxiliary 
line. The minimum of T takes place in the interval 
lgT c „ w 8.5-8.7, because at 8.0 < lgT c „ < 8.0 + 0.5 
the neutrino energy losses are partly compensated by 
the latent heat release. The lowest temperatures 
are realized at T cp 3> 10 s K, when the nucleon heat 
capacity suffers the strongest suppression. 

On vertical parts of isotherms below the lower dashed 
line the cooling is regulated by the neutron super- 
fluidity alone. With growing T cn in this domain the 
temperature T{ varies in a different manner than with 
growing T cp in the domain of the purely proton su- 
perfluidity. This happens because at lgT c „ > 7.9 
the neutrino emission due to the neutron pairing is 
important. It speeds up NS cooling in the range 
7.9 < lgT cn < 7.9 + 0.7 (see above). If 7.9 < lgT c „ < 
7.9 + 0.5 the neutrino cooling is partly compensated by 
the latent heat release. Therefore the minimum of T$ 
takes place in the interval 7.9 + 0.5 < lgT c „ < 7.9 + 0.7. 
It is not so deep as in the upper part of Fig. ^ since 
the nucleon heat capacity is suppressed only partly. 

A strong neutron superfluidity reduces the heat ca- 
pacity stronger and the neutrino luminosity weaker 
than a strong proton superfluidity. Owing to the weak- 
ness of the neutrino energy losses this difference is suf- 
ficient for a NS with high T c „ and normal p to cool 
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Table 2: Observational data 



Source 


ig* [yr] 


Model V 


lgT- [K] 


Confidence 
level a > 


Reference 


RX J0822-43 


3.57 


Hydrogen atm. 


a 90+O.O2 
°- zo -0.02 


95.5% 


Zavlin et al. (1999b) 


Black body 


6 fil+ 005 


95.5% 


Zavlin et al. (1999b) 


IE 1207-52 


3.85 


Hydrogen atm. 


6 10+ 005 


90% 


Zavlin et al. (1998) 


Black body 


a 40+0 02 

D - 4y -0.01 


90% 


Zavlin et al. (1998) 


RXJ0002+62 


3.95 c > 


Hydrogen atm. 




95.5% 


Zavlin et al. (1999a) 


Black body 


6.i8lg:iS 


95.5% 


Zavlin et al. (1999a) 


PSR 0833-45 
{ Vela ) 


4 4 d ) 


Hydrogen atm. 


5.90±85* 


90% 


Page et al. (1996) 


Black body 


R 9 4 +0.03 
D - z ^-0.03 




Ogelman (1995) 


PSR 0656+14 


5.00 


Hydrogen atm. 


r- 7 9+0.04 
°- ,z -0.02 




Anderson et al. (1993) 


Black body 




90% 


Possenti et al. (1996) 


PSR 0630+178 
(Geminga) 


5.53 


Hydrogen atm. 


r 9 r+0.08 

°- ZJ -o.oi 


90% 


Meyer et al. (1994) 


Black body 


r 7C+0.05 
°- ,o -0.08 


90% 


Halpern, Wang (1997) 


PSR 1055-52 


5.73 


Black body 






Ogelman (1995) 



a ' Confidence level of (90% and 95.5% correspond to 1.64c and 2<r, respectively); dash means that the level is not indicated 
in cited references. 

Model used for interpretation of observation. 
c ) The mean age taken according to Craig et al. (1997). 
^' According to Lyne et al. (1996). 



in a different way than at equally high T cp and nor- 
mal n. Strong neutron superfluidity delays the cooling 
of those NSs which would be at the neutrino cooling- 
stage or at the neutrino-photon transition stage in case 
they were nonsuperfluid. This is demonstrated in Fig. 
H (for PSR 0656 + 14); in the absence of superfluidity, 
at t — 10 5 yr this pulsar would be at the transition 
stage. Strong neutron superfluidity accelerates cooling 
of older NSs, e.g., of the Geminga's age (Fig. ||). 

Vertical parts of isotherms in Fig. || intersect the 
auxiliary line at temperatures T cp , at which the proton 
superfluidity appears. This superfluidity leads to the 
latent heat release and to the growth of Ti in a strip 
of width ~0.5 in lgT cp above the lower dashed line. 
With further growth of T cp the heat capacity is strongly 
reduced and cooling is accelerated. For a very strong 
joint superfluidity of n and p (T cn 3> 10 7,9 K, T cp ;§> 
10 7 ' 9 K, the very right upper corner of Fig. ||), the 
nucleon heat capacity and the neutrino luminosity of 
the NS core are fully suppressed, and the cooling is 
governed by the electron heat capacity. 

Figure || is analogous to Fig. g, but corresponds to 



a younger NS. Its neutrino luminosity is somewhat 
higher and the relative contribution of neutron-pairing 
neutrinos is smaller. Neutrino emission produced by 
pairing affects the cooling weaker. 

3 Comparison of calculations and 
observations 

3.1 Observations and their interpreta- 
tion 

By now, thermal X-ray radiation has been detected 
from seven isolated NSs (Table ||). All the sources 
are reliably identified. Four of them (Vela, Geminga, 
PSR 0656+14 and PSR 1055-52) are radiopulsars. The 
other three NSs (IE 1207-52, RX J0002+62 and RX 
J0822-43) are radio silent, but are observed as point- 
like sources of soft X-rays in supernova remnants. 
Their identification as isolated NSs is confirmed by ob- 
servations of coherent pulsations of X-ray emission with 
NS-spin periods (excluding IE 1207-52) and smooth 
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lgt[yr] 

Figure 7: Observational data on surface tempera- 
tures of NSs obtained (Table 2) in the models of black 
body radiation (filled circles) and hydrogen atmosphere 
(open circles). Shaded regions show ranges of T s °° filled 
by standard (diagonal shading) and enhanced (horizon- 
tal shading) cooling curves of NSs with different T cn 
and T cp (from 10 6 to 10 10 K). The solid line shows 
standard cooling (1.30 M Q ) at specified values lgT cra 
and lg T cp (given near the curve). Dotted lines ex- 
hibit standard (upper line) and enhanced (lower line) 
cooling of a non-superfluid NS. The dashed and dot- 
and-dashed lines show, respectively, standard and en- 
hanced cooling of a NS which possesses an envelope of 
mass 7-10~ 10 Mq composed of light elements (see Sect. 
3.3). 

(modulation depth < 15 — 40%) light curves, and also 
by small ratios of the optical to X-ray luminosities, 
which are not typical for accreting NSs in binaries. 

As a rule, an X-ray spectrum from an isolated NS 
contains two components. The spectrum of the softer 
component is similar to the black body spectrum for a 
temperature T s w 10 5 -10 6 K. According to the cooling 
theory, the surface of a middle-aged NS, t w 10 3 -10 6 
yr, should have about the same temperature. It is as- 
sumed, therefore, that the main contribution into the 
soft component comes from the thermal surface emis- 
sion. The hard component is associated either with 
a thermal emission from hot polar caps on the pulsar 
surface (heated by an inverse current of energetic par- 
ticles from the magnitosphere) or with a non-thermal 



magnetospheric emission. 

The thermal NS radiation is often fitted by the black 
body spectrum. However, the interpretations based on 
the spectra, which are obtained by modelling of the NS 
atmospheres, are more realistic. At the same effective 
temperature T s the black body spectrum appears to be 
softer. Therefore, while interpreting the same observa- 
tions, the "black body" temperature T s is noticeably 
(typically, by a factor about 1.5-3) higher than the 
"atmospheric" temperature. 

Construction of the atmosphere models in a prac- 
tically important range of temperatures and surface 
magnetic fields B has not been completed yet. Only 
the hydrogen-helium and iron atmosphere models with 
weak magnetic fields B < 10 9 G are really reliable as 
well as the hydrogen-helium models with strong fields 
(B > 10 11 G) for T s > 10 6 K. If these conditions are 
fulfilled, the thermal radiation spectrum and the NS 
light curves are better described by the atmospheric 
model than by the black body model. Physical param- 
eters inferred from atmospheric interpretations (radi- 
ating surface area, distance to a NS, column density of 
interstellar gas, etc.) are in better agreement with the 
data obtained from independent observations in other 
wavelengths (see, e.g., Page et al. 1996, Zavlin et al. 
1998, 1999a, b). The black body interpretation often 
gives less realistic values of the parameters and meets 
some difficulties in explaining all sets of observational 
data. 

In Table | the NSs, the sources of the thermal emis- 
sion, are ordered according to their ages. For every NS, 
we present its age, the method of interpretation of ther- 
mal radiation, and the effective surface temperatures 
T s °° = T s y/1 - Rg/R obtained by a given method. For 
the pulsars, we present the dynamic ages, while for the 
radio silent objects we give the ages determined from 
the associated supernova remnants. Although the val- 
ues of M and R (and, therefore, of the gravitational 
redshift) used in our models (Sect. 2) differ slightly 
from those obtained (or adopted) from interpretation 
of observations, the temperatures T s °° from Table || can 
be compared with our cooling calculations due to the 
following reasons. First, our calculations are not very 
sensitive to variations of NS mass or radius (see Sect. 
3.3 below). Second, the confidence ranges of M and/or 
i?, obtained (adopted) from interpretation of observa- 
tions, include, as a rule, our theoretical values. 

Among the atmospheric interpretations presented in 
Table [| the results for IE 1207-52, RX J0002+62, 
RX J0822-43 and for the Vela pulsar seem to be more 
reliable. The advantage of interpretation of their spec- 
tra by the hydrogen-helium atmosphere models is ad- 
ditionally confirmed by the fact that fitting of ob- 
servations by the iron atmosphere models or by the 
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Figure 8: Confidence regions of T cn and T cp which correspond to the "black body" surface temperatures (Table 
2, Fig. |) of RX J0002+62, PSR 0656+14 and Geminga in the models of enhanced (left) and standard (right) 
NS cooling. 



black body spectrum is of lower statistical significance. 
The Geminga and PSR 0656+14 pulsars are older and 
cooler. The important contribution into the spec- 
tral opacity of their atmospheres should come from 
the effects of motion of neutral and partly ionized 
atoms across magnetic fields (Potekhin & Pavlov 1997, 
Bezchastnov et al. 1998). However, only approximate 
atmosphere models have been used for interpretation so 
far, in which these effects have been ignored. The pres- 
ence of the hydrogen atmosphere of the Geminga pulsar 
is additionally confirmed by a possible discovery of the 
proton cyclotron line in the optical spectrum (Bignami 
et al. 1996, Martin et al. 1998). The "atmospheric" 
temperatures of the Geminga and PSR 0656+14 pul- 
sars are likely to be closer to the reality than the black 
body temperatures, but less reliable than for younger 
NSs. 

We will assume that the atmospheric and black body 
interpretations are equally possible and analyze the 
values of T cn and T cp which can be inferred from these 
data. 

3.2 Necessity of superfluidity for inter- 
pretation of observations 

In Fig. the cooling curves are compared with the 
data of Table 0. Diagonal and horizontal shadings show 



the ranges of the surface temperature T£°(i) filled, re- 
spectively, by different standard and enhanced cooling 
curves calculated by varying T cn and T cp from 10 6 to 
10 10 K. Dashes show cooling of nonsuperfluid stars. 

The "nonsuperfluid" curves are seen to be in poor 
agreement with observations. On the other hand, the 
observations can be explained by assuming superflu- 
idity in the NS cores. This is illustrated (Fig. 0) by 
the standard cooling curve (the solid line). The val- 
ues of T cn and T cp are chosen in such a way to hit the 
maximum number of observational points at once. 

According to Fig. 0, all the "atmospheric" temper- 
atures as well as the "black body" temperatures of 
RX J0002+62, PSR 0656+14 and Geminga hit the al- 
lowed regions of the standard and enhanced cooling of 
superfluid NSs. Thus our cooling calculations can be 
compared with the "atmospheric" and "black body" 
temperatures T s of these sources. 

High black body surface temperatures of RX J0822- 
43, IE 1207-52 (not presented in Fig. but given m 
Table 0), Vela, and PSR 1055-52 are not explained 
by our models but can be explained by other mod- 
els of cooling NSs, including the models with super- 
fluid cores. For instance, observations of the Vela pul- 
sar agree with the standard cooling of a superfluid NS 
(T cn = 10 7 K,T cp = 10 10 K) possessing an outer en- 
velope of mass ~ 10~ 9 M Q composed of light elements 



10 




Ig T cn [K] lg T cn [K] 

Figure 9: Same as in Fig. but for the enhanced cooling and "atmospheric" interpretation (Table 2) of spectra 
of thermal radiation from RX J0822-43, IE 1207-52, RX J0002+62, Vela, PSR 0656+14 and Geminga. On the 
right panel we show in more detail the region (shaded rectangle) in which the allowed values of T cn and T cp of 
five latter objects are either close or intersect. 



(Potekhin et al. 1997). High black body temperatures 
of RX J0822-43, IE 1207-52 and PSR 1055-52 may in- 
dicate the presence of some additional reheating mech- 
anism inside these sources (see, e.g., Umeda et al. 1993) 
or the presence of superstrong (B > 10 14 G) magnetic 
fields (Heyl & Hernquist 1997). Finally, one cannot ex- 
clude the possibility that the black body interpretation 
of their spectra is incorrect (see above). 

3.3 Results 

As seen from Fig. Q by assuming the presence of nu- 
clcon superfluidity in the NS cores, we can explain the 
majority of NS observations either by standard or by 
enhanced NS cooling. These observations include six 
"atmospheric" interpretations (RX J0822-43, IE 1207- 
52, RX J0002+62, Vela, PSR 0656+14, Geminga) and 
three "black body" ones (RX J0002+62, PSR 0656+14, 
Geminga). For simplicity, we assume that the inter- 
nal structure of these NSs is the same, and, particu- 
larly, the NSs have the same mass. Then the critical 
temperatures of nucleons in their cores are the same. 
Let us analyse the confidence regions of T cn and T cp 
constrained by observations. Including observations of 
several NSs allows us to reduce these regions. 

The regions in question are plotted in Figs. p| |l0|. 



Figure |8J corresponds to the standard and enhanced 
cooling of NSs (with masses 1.30 M and 1.48 M , 
respectively) with the black body spectrum. Figure 

corresponds to the enhanced cooling of the stars 
(1.48 Mq) possessing hydrogen atmospheres. Finally, 
the standard cooling of NSs (1.30M Q ) with hydrogen 
atmospheres is shown in Fig. [l0| In every figure, the 
lines of different type enclose the confidence regions of 
T cn and T cp associated with the error bars of the ob- 
served NS surface temperatures T s °° (Table 2). Corre- 
spondence of the lines to the selected NSs is displayed 
in Fig. g, For the PSR 0656+14 and Geminga pulsars, 
isotherms are taken from Figs. In each figure the 
actual overall confidence region of T cn and T cp lies at 
the overlap of the confidence regions of all objects. 

Figures |s|-|Io| show that observations of several NSs 
at once can be explained either by the standard or by 
the enhanced cooling, adopting either black body or 
atmospheric interpretations of the spectra. In all the 
cases there are the ranges of T cn and T cp close or joint 
for all NSs; they do not contradict to the microscopic 
theories of nucleon superfluidity in NSs (see Takatsuka 
& Tamagaki 1997 and references therein). 

According to the left panel of Fig. |[ by adopting 
the enhanced cooling and the black body interpreta- 
tion of observations we obtain two confidence regions 
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of T cn and T cp ; each explains cooling of three objects 
at once. The first region corresponds to a moder- 
ate neutron superfluidity (lgT cn « 8.1) and a strong 
proton superfluidity (lg T C p w 9.5); the second, wider 
region corresponds to a strong neutron superfluidity 
(lgT c „ « 9.75-10.0) and a moderately weak proton su- 
perfluidity (lg T cp < 8.3). For the standard cooling of 
a NS with the black body spectrum (the right panel 
of Fig. ||) there are also two regions of T c „ and T cp for 
the three NSs, but they are somewhat different than for 
the enhanced cooling. The first region corresponds to a 
moderate neutron superfluidity (lgT c „ = 8.0 8.1) and 
a moderately strong proton superfluidity (lgT cp = 8.2- 
9.0), while the second one is associated with a strong 
neutron superfluidity (lgT cn > 9.8) and a moderately 
weak proton one (lgT cp < 8.2). 

In the case of the enhanced cooling and the "at- 
mospheric" surface temperatures (Fig. ||) there is the 
only region of joint or very close values of T cn and 
T cp for five NSs (the shaded rectangle on the right 
panel). It corresponds to a moderate neutron super- 
fluidity (lgT c „ ks 8.2-8.35) and a strong proton su- 
perfluidity (lgT cp ~ 9.45-9.65). This superfluidity 
(lgT cp = 8.3, lgT C p = 9.6) allows us to describe the 
sixth object from the "atmospheric" set (RX J0822- 
43), if we assume that it possesses a thermally insulat- 
ing envelope of light elements with the mass equal to 
7 • 10~ 10 M Q (see Fig. ^). However even in the absence 
of the envelope the confidence region of T cn and T cp for 



this object is sufficiently close to the joint confidence 
region indicated above. 

Finally, for the standard cooling and the "atmo- 
spheric" spectrum (Fig. [!(]), there is again the only 
confidence region of T c „ and T cp , where the critical 
temperatures are nearly the same for the five NSs; it 
corresponds to a moderately strong superfluidity of n 
and p (lgT cn w 8.25-8.35, lgT cp « 8.0-8.4, shaded 
rectangle). The observations of RX J0822-43 can also 
be explained by the presence of the same superflu- 
idity (lgT c „ = 8.3, lgT C p = 8.2) assuming that the 
object possesses an envelope of light elements of the 
7 ■ 1O _1O M mass (see Fig. 0). In the absence of the 
envelope, the confidence region of T cn and T cp lies out- 
side of the joint confidence region since the error bar of 
T s °° for this object is situated at the boundary of the 
shaded region (filled by the cooling curves of supernuid 
NSs). 

Therefore, the constraints on the critical tempera- 
tures T c „ and T cp depend on cooling type and on inter- 
pretation of the NS thermal radiation. The fact that 
the NSs have actually different masses has little ef- 
fect on our results as long as we assume that T cn and 
T cp are independent of density. The standard cooling- 
curves of NSs with masses lower than the threshold 
mass M cr = 1.442 Mq (of switching on the direct Urea 
process) practically coincide. The same is true for the 
enhanced cooling curves of NSs (M > M cr ) as long as 
their cores contain a moderate or strong superfluidity 
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in the NS cores, one needs new cooling simulations (for 
different model equations of state and nucleon critical 
temperatures taking into account variation of T cn and 
T cp over the NS cores), new high-quality observations 
of NS thermal emission and new atmospheric models 
for interpretation of observations. 

We arc grateful to G.G. Pavlov and V.E. Zavliv 
for providing us with the data on RX J0822-43 and 
RX J0002+62 prior to publication. The work is partly 
supported by the grants of RFBR (96-02-16870a), DFG- 
RFBR (96-02-00177G) and INTAS (96-0542). 
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of n and/or p with T c > 10 8 K. With weakening of 
the superfluidity the difference between enhanced and 
standard cooling curves increases reaching maximum 
for nonsupcrfluid NSs. 

The existence of the same critical temperatures for 
several NSs at once is quite unexpected. We have ex- 
pected different joint confidence regions of T cn and T cp 
for different pairs of NSs. The result is even more 
surprising taken into account simplicity of our cool- 
ing models. It would be interesting to confirm (or re- 
ject) these results using advanced cooling models and 
a larger number of objects. 

As a whole, one needs stronger superfluidity of at 
least one nucleon species for the enhanced cooling mod- 
els. In the case of the black body interpretation of NS 
radiation, the stellar cores may contain a strong neu- 
tron superfluidity, T cn « 10 9 7 -10 10 K. In neither case 
one needs strong superfluidity of n and p simultane- 
ously to explain the observations. This indicates that 
the equation of state of the NS cores cannot be too soft 
(the softness would mean weak nucleon-nucleon repul- 
sion at small separations which would induce especially 
strong pairing). 

Finally let us mention that we can satisfy observa- 
tions by varying the only parameter, T cp , if we as- 
sume the presence of neutron superfluidity with T cn « 
10 81 -10 83 K. For the standard cooling, this parameter 
should lie in the moderate range 10 8 -10 9 K, while 

for the enhanced cooling it should be larger, 10 9 45 - 
10 9.65 K 

4 Conclusions 

Our results show that cooling of NSs is very sensitive to 
superfluid state of neutrons and protons in their cores. 
This confirms the possibility (Page & Applegate 1992) 
to explore the nucleon superfluidity in the NS cores, 
the fundamental property of matter of supranuclear 
density, from NS cooling. Our calculations confirm also 
that the neutrino emission due to nucleon pairing plays 
important role in NS cooling. 

We have used the simplified NS models with two 
fixed masses, 1.48 M and 1.30 Mq, and nucleon criti- 
cal temperatures constant over the stellar cores. Even 
in this case we have been able to explain the major- 
ity of observations of the thermal radiation from Vela, 
Geminga, PSR 0656+14, RX J0002+62, IE 1207-52 
by four models: by the standard or enhanced cooling 
of NSs using the black body or "atmospheric" inter- 
pretations of their spectra. For every model, we have 
obtained its own, sufficiently constrained confidence re- 
gions of the critical temperatures T cn and T cp . 

Our results are preliminary. For obtaining more reli- 
able conclusions on superfluidity and equation of state 
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